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Research on graphene has revealed remarkable phenomena arising in the honeycomb lattice. However, 
the qnantum spin Hall effect predicted at the K point could not be observed in graphene and other 
honeycomb structures of light elements dne to an insufficiently strong spin-orbit coupling. Here we show 
theoretically that 2D honeycomb lattices of HgTe can combine the effects of the honeycomb geometry and 
strong spin-orbit coupling. The condnction bands, experimentally accessible via doping, can be described 
by a tight-binding lattice model as in graphene, but including mniti-orbital degrees of freedom and spin- 
orbit coupling. This resnlts in very large topological gaps (np to 35 meV) and a flattened band detached 
from the others. Owing to this flat band and the sizable Coulomb Interaction, honeycomb structures of 
HgTe constitute a promising platform for the observation of a fractional Chern insniator or a fractional 
quantum spin Hall phase. 
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The discovery of graphene has confronted us with a ma¬ 
terial that exhibits fascinating electronic propertieJ^I such as 
zero-mass carriers, persisting conductivity at vanishing den¬ 
sity at the Dirac poinP, Klein tunnellinj^ and an anomalous 
quantum Hall effecP®. Nevertheless, the absence of a band 
gap in its spectrum prevents its use as a field-effect transistor, 
and its weak spin-orbit coupling (SOC) hampers the possibil¬ 
ity to realize the quantum spin Hall effect (QSHEj^ and use it 
for quantum spintronics. The prospect of artificial graphene 
samplers that display the lacking properties has motivated 
research in various types of honeycomb lattice, such as ar¬ 
rays of ultracold atoms^ rnolecular graphen^i^, organometal- 
lic latticeJ ^, and 2D electron gases subject to a geometric ar¬ 
ray of gateJ®3. Recently, an alternative path came from the 
self-assembly of semiconductor nanocrystals forming atom¬ 
ically coherent 2D structures with a long-range honeycomb 
pattern; the thickness and honeycomb period are defin ed by 
the size of the nanocrystals, which is in the range of 5 nrrP32l 
Honeycomb lattices of PbSe and CdSe nanocrystals have been 
fabricated, with astonishing atomic coherence due to the ori¬ 
ented attachment of the nanocrystals. Theoretical investiga¬ 
tions have shown that CdSe superlattices formed in such a 
way exhibit Dirac cones at two energies and nearly disper¬ 
sionless bands. Unfortunately, these flat bands are connected 
to the nearby higher energy bands, and the SOC gaps in the 
conduction band are very smalP^ 

In this work, we propose a design for robust topological in¬ 
sulators that combine the properties of the honeycomb lattice 
and strong SOC. We consider three different types of HgTe 
layers with superimposed honeycomb geometry and present 
atomistic tight-binding (TB) calculations of their conduction 
band structure that we accurately describe by a 16-band effec¬ 


tive model. Such lattices take advantage of the multi-orbital 
degrees of freedom in the honeycomb setup, allied to the 
strong SOCP. The access to multi-orbital degrees of freedom 
allows for further manipulation of the topological properties. 
While the Haldan^SI and Kane-Mel^ models originally con¬ 
cerned honeycomb lattices characterized by a single orbital 
per site (e.g., in graphene) and isotropic nearest-neighbour 
(NN) hopping integrals between them, multi-orbital models 
have attracted much attention recently, in particular because 
they seem to be a paradigm to generate topologically non¬ 
trivial flat bandiiSJRlI. The flat band structure opens the way 
to the realization of interesting strongly correlated states, such 
as fractional QSHE, fractional Chern in sulators, or ferromag¬ 
netic fractional Chern insulatorj22lIl2lllII] xhe currently pro¬ 
posed HgTe lattices exhibit conduction bands characterized 
by large topological gaps and an isolated flat band. We con¬ 
clude that, depending on the position of the Eermi level, not 
only QSHE could be observed in these structures, but also 
fractional QSHE or fractional Chern insulator phases, as on¬ 
site and NN Coulomb-interaction parameters are found in the 
energy range required for their realization. 


Results 

Design of the HgTe honeycomb lattices. The three struc¬ 
tures that we consider have in common that the [111] direction 
of the zincblende lattice is perpendicular to the plane. The first 
system is ins pired by recent results obtained for nanocrystal 
self-assemblySEIl j e., a graphene-type superlattice of trun¬ 
cated cubic nanocrystals attached via (110) facets. The sec¬ 
ond system consists of spheres connected by short cylinders, 
which allow us to vary the coupling between honeycomb lat- 
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FIG. 1. Nanocrystal lattices and their conduction band dispersions, a, honeycomb nanoribbon formed by the HgTe (CdSe) nanocrystals. 
Hg (Cd) atoms are in yellow, Te (Se) atoms are in grey. Each nanocrystal has a truncated nanocube shape, the vertices of which are given by 
all permutations of [±1, ±(1 — q), ±(1 — q)]l, where q is the truncation factor and 21 is the size of the original nanocube before truncation. 
The honeycomb lattice spacing a, i.e., the center-to-center distance between neighbor nanocrystals, is defined as a = (2N + l)ao/v/2 where 
N is an integer and ao is the cubic lattice parameter of HgTe (CdSe). The nanocrystals are attached via (110) facets (a = \/2(2 — q)l). The 
arrows along the ribbon indicate the electron propagation in the helical edge states present in the quantum spin Hall phase. Red and blue colors 
correspond to top and bottom edge for spin up, bottom and top edge for spin down, respectively. b,c, band dispersions for the bulk resulting 
from the atomistic tight-binding calculations with q = 0.5 (a = 5.0 nm for HgTe, 4.7 nm for CdSe). d, same for the HgTe superlattice 


resulting from the effective model. Insets in (c,d) show the s bands in the gap region with 


higher magnification ^2K/3 = ^^)]- 


tice sites in a convenient way. The third system corresponds 
to a honeycomb array of cylinders. These two last systems 
may be experimentally realized by gas-phase deposition and 
lithography. We will show below that the multi-orbital topo¬ 
logical effects are common to these three different structures. 

The HgTe superlattices proposed here differ fundamentally 
from HgTe/CdTe quantum wells, where the QSHE has been 
predicted^ and experimentally observedSI In the latter sys¬ 
tem, the appearance of the QSHE at the T point is connected 
to band inversion and vanishes for quantum wells of thickness 
below 6 nm. Here, instead, the effect occurs at the K point, 
and it is driven by the honeycomb nanogeometry, allied to the 
strong SOC of the composing HgTe nanocrystals. This dis¬ 
tinction is important because zero modes (Majoranas) bind to 
topological lattice defects w hen th e band gap opens at a non- 
r point in the Brillouin zondSEl and hence, the underlying 
topological order can b e dete cted by measuring the structure 
of the topological defecP^I^Il. 

Band structure of lattices of HgTe nanocrystals. To un¬ 
veil the topological properties of these systems, we have per¬ 
formed atomistic tight-binding (TB) band-structure calcula¬ 
tions. We use a basis of twenty atomic orbitals on each atom 
of the nanocrystals, including the spin degree of freedom. The 
methodology is described in Ref. and is summarized in 
the Methods section. In comparison to bulk HgTe, the elec¬ 
tronic structure of HgTe superlattices is characterized by a 
large bandgap due to the strong quantum confinement. Con¬ 
duction and valence bands close to this gap are composed of 
several minibands and minigaps due to the periodic scattering 
of the electronic waves in the honeycomb structurJiSI in the 


following, we only discuss the physics of the sixteen lowest 
conduction bands (Fig.[2l. 

A typical honeycomb lattice composed of HgTe nanocrys¬ 
tals is shown in Fig.[^ and the related conduction-band struc¬ 
ture is displayed in Fig. [T};. In order to visualize the effects of 
SOC, we compare it to the band structure of the same honey¬ 
comb structure composed of CdSe (Fig.[^). The strong SOC 
in HgTe gives rise to effects which are absent in CdSe^^ The 
overall behaviour of the band structure can be understood as 
follows. Each individual nanocrystal is characterized by two 
states with s envelope wave-function and six p states at higher 
energy. In the honeycomb structure, strong coupling between 
the wave-functions of neighbour nanocrystals leads to the for¬ 
mation of sixteen bands grouped into two manifolds of four 
(s) and twelve (p) bands that are well separated. The s bands 
have the same type of dispersion as the tt and tt* bands in real 
graphenJil In the case of CdSe (Fig.[^), these bands are spin- 
degenerate and are connected at the K and K' points of the 
Brillouin zone, where their dispersion is linear (Dirac points). 
In HgTe, instead, the s bands exhibit a small gap (5.7 meV) 
at the K points and have a quadratic dispersion (see Fig. 0)- 
In addition, they are characterized by a visible spin splitting 
at all points of the Brillouin zone, except at T and M. Among 
the twelve p bands of CdSe nanocrystal superlattices, eight 
have a small dispersion and the other four basically behave 
like the (Dirac) s bands. Four flat bands are built from the p^ 
states perpendicular to the lattice, which are not very disper¬ 
sive because pz-pz (tt) interactions are weak. Four other bands 
(Px,y), respectively above and below the p Dirac band, are flat 
due to destructive interferences of electron hopping induced 
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by the honeycomb geomettyE^^. In honeycomb lattices of 
HgTe nanocrystals, the SOC induces spin splitting, opens a 
large gap at K in the p-like Dirac bands, and produces a con¬ 
siderable detachment of the lowest flat p-band from the Dirac 
p-band (Fig.[T]:). The effects of the SOC are so strong that it 
is hardly possible to recognize the Dirac bands. 

These unexpected features allow for the realization of sev¬ 
eral topological states of matter, by doping the system using 
a field-effect transistor or electrolyte gel gatinj^. At zero en¬ 
ergy, the undoped system is a semiconductor with a trivial 
gap of about 0.4 eV. Upon doping the material with one elec¬ 
tron per nanocrystal, the small s-like topological gap may be 
reached, whereas for fillings between 2 and 3 electrons per 
nanocrystal, the fractional quantum (spin) Hall regime may 
be realized at the flat band. For a doping level of 4 elec¬ 
trons/nanocrystal one reaches the QSHE gap. At this point, 
we should emphasize that doping of nanocrystals with up to 
10 electrons has already been demonstrated experimentalljES 
therefore, all the interesting regimes that we discuss are at 
reach with the existing technology, at a simple switch of the 
doping level. Other examples of band structures for lattices 
of HgTe nanocrystals with different size or truncation factors 
are presented in Supplementary Fig. 1 and are discussed in 
Supplementary Note 1. They all show large topological gaps, 
especially in the p sector. 

The topological properties of the bands are most trans¬ 
parently studied through an edge-state analysis in a ID 
nanoribboiPE. We consider a zigzag ribbon composed of six¬ 
teen nanocrystals (34740 atoms) per unit cell. Figure]^ shows 
that edge states are crossing the three gaps between the p 
bands as well as the gap between the s bands. These results 
also hold for armchair ribbons. In Supplementary Fig. 2a, we 
present the band structure for another nanoribbon, which has 
two inequivalent edges. Still, helical edge states characteris¬ 
tic of the QSHE are found, as shown by the 2D plots of the 
wavefunctions (Supplementary Eigs. 2b-e). 

Effective model. The band structures resulting from 
the atomistic TB calculations are accurately described by a 
sixteen-band effective model (Eig. where each nanocrys¬ 
tal is treated as one site on a honeycomb lattice. The effec¬ 
tive TB model is written in the basis of the four aforemen¬ 
tioned orbitals {s,Px,Py,Pz) per site as Tfeff = + iTNN + 

^^ISO+^^RSO- Here, incorporates the on-site energies Eg, 
Ep^ = Ep ^, and Ep^. The edge nanocrystals have a slightly 
different value of Eg compared to the bulk, in order to account 
for the different number of neighbours. The term 

,a, ( 1 ) 


encodes NN hopping, where {i,j) denotes NN sites, a = f, i 
denotes spin, and b and b' the orbitals. The coupling coeffi¬ 
cients Vi^b-,j,b' are expressed in terms of the hopping param¬ 
eters Vgsa, Vppa, VppTT, and Vgpa, following the notations of 
Ref.gll 

The intrinsic SOC term Tfiso couples the electron orbital 
angular momentum L and spin S = <t/ 2. In the p sector, it is 
encoded through the on-site term Af 3 QL • cr. There is no on¬ 


site term in the s sector because the orbital angular momen¬ 
tum is “frozen”. Eor the same reason, in graphene, the on¬ 
site intrinsic SOC term is absent because the hybridiza¬ 
tion freezes the orbital momentum in the p^ state. However, 
as shown by Kane and Mel^, the intrinsic SOC introduces a 
next-nearest-neighbour (NNN) hopping term, which is written 
as 

^ISO = iAjso ^ ^ ^ (2) 

Here, the summation is over NNNs, and = ±1, with the 
sign depending on the outer product of the two NN vectors 
that connect sites i and j. The Rashba SOC term, proportional 
to the cross product of momentum and spin, p x S, is written 
as a NN-hopping term 

^^RSO = i EE ^ Uj )]a/3Cj,b',p. (3) 

(ij) 

b,b' 

The coupling coefficients Ji^b-JM' have the same structure as 
the Vi^b-j,b' for the ordinary NN hopping and are expressed in 
terms of jgga, Jppa, and The sp term may be neglected 
due to the large energy separation between the s and p bands. 

In Table [Ij we present typical values for the parameters 
obtained numerically using least-squares fitting to the band 
structure of Eig. The band structure of the effective model 
is shown in EigTm. As expected, the on-site term Ajqq is 
much larger than me NNN term AfgQ, explaining the open¬ 
ing of a very large gap at the K point and the detachment of 
the flat bands in the p sector. Using the effective model, the 
non-trivial topology of the bands is confirmed by the calcu¬ 
lation of the Z 2 topological invariant, the spin Chem number 
(Methods). 

The effective model yields a band structure for the ribbon 
(Fig . 1^) in excellent agreement with the atomistic TB calcu¬ 
lations (Eig.[^). The red and blue colours in Eigs.|^,b indi¬ 
cate the expectation value {yaz) = where y is 

the coordinate perpendicular to the ribbon edges. This expec¬ 
tation value allows us to identify helical edge states that come 
in pairs with identical dispersion, opposite spin, and live on 
opposite edges. Both atomistic and effective TB models show 
that the gaps in the s and p sectors exhibit helical edge states, 
characteristic of the QSHE. 

When the Rashba coupling is neglected in the effective 
model, all states are spin degenerate. The Rashba term in¬ 
duces a small splitting in energy, and tilts the spins slightly 
away from the perpendicular direction. In Eigs. |^, we plot 
the spin direction on each site of the zigzag ribbon for four 
sets of four edge states (two in each edge). The spin direction 
is always perpendicular to the edge, i.e., the spin lies in the yz 
plane if we choose the x direction to be parallel to the edge. 
The localization of the selected states on the edges is visible 
from the weight of the wave function, indicated by the size of 
the circles in the figure. The colours of the circles are deter¬ 
mined by the local value of yaz, and correspond to the colours 
in the dispersion (Eig.|^,b). The site dependence of the spin 
direction leads to interesting spin textures. Eor the edge states. 
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FIG. 2. Topological edge states and non-trivial gaps in honeycomb lattices of HgTe nanocrystals, a, conduction bands calculated using 
the atomistic tight-binding method for a zigzag ribbon composed of sixteen nanocrystals per unit cell {q = 0.5, body diagonal of 5.0 nm, cell 
length I = 8.7 nm). b, same but computed from the effective Hamiltonian, c, spin orientation on each site for a selection of states calculated at 
k = 3l7r/(301). A vertical arrow indicates that the spin is along the z direction, perpendicular to the lattice. The size of the circles represents 
the weight of the wave function on each site. In each figure, the colour indicates the expectation value {yo^), i.e., red and blue correspond 
to top and bottom edge for spin up (bottom and top edge for spin down), respectively. At each energy E, there are two states which live on 
opposite edges with opposite spin, therefore they are represented by the same colour. The bulk states are grey. 


on-site NN hopping Rashba SOC intrinsic SOC 

^buik _ Q 305 gV Vssa = —26.4meV ysso- = 0.56 meV Afso = 0.71 meV 

^edgo _ Q 37QeV Vpptr = 45.6 meV jppa = 1.50 meV AjgQ = 15.8 meV 

Ep^ = 0.691 eV Vpp-n = —2.7meV 7 pp,r = 0.80 meV 

Ep^ = 0.691 eV Vapa = 31.1 meV 
Ep^ = 0.747eV 


TABLE I. Parameters of the effective model. Parameters derived for the lattice of HgTe nanocrystals described in Fig.fT] Eg, Ep ^, Ep ^, and 
Ep^ are the on-site energies on the s, Px, Py, and p^ orbitals, respectively. In ribbons, the edge nanocrystals have a slightly different value of 
Eg compared to the bulk (75^“*’' 7 ^ Vgga, Vppa, Vppn, and Vgpo- are the hopping parameters, following the notations of Ref. [42 Tsstr, 

7 ppcr, and 'ypp-w are the terms describing the Rashba SOC, following the same notations. The intrinsic SOC is defined by Afso and AfgQ on s 
and p orbitals, respectively. 


the typical spin texture is almost smooth: Going from one 
edge to the other, the spin direction interpolates between (al¬ 
most) up and (almost) down in a rotational manner. A slight 
tilt of ^ 3° is observed for the edge states in the p bands. In 
the s bands, the tilt is stronger, similarly to graphenJS] Here, 
the spin vector at the edge site points ~ 30° away from the 
vertical. This difference in tilt can be explained by the much 
larger intrinsic SOC in the p sector than in the s one, whereas 
the two Rashba couplings jssa and 'yppa are of similar order 
of magnitude. Strictly speaking, one should denote this state 
a Z 2 topological insulator, but since the Rashba SOC is ex¬ 
tremely small, one can think of an approximate QSHE. 

Flat band and Coulomb interactions. The large gaps in 
the p sector are mainly due to the intrinsic SOC, which, con- 
trarily to the Kane-Mele model, is described by an on-site 
term (AfgQ = 15.8 meV). In particular, the intrinsic SOC sep¬ 
arates the lowest flat band from the other p bands, with large 
gaps, e.g., 35 meV in the case of Fig. [^. This gap ranges 
from 13 meV to 36 meV when we vary the nanocrystal size 


S^out 

U V 

Band width 

Absolute gap Gap at F 

14 14 

14 6 

48 meV 23 meV 

76 meV 43 meV 

20 meV 

35 meV 56 meV 


TABLE 11. Coulomb energies. On-site (U) and nearest-neighbor 
(C) Coulomb interaction energies calculated for a honeycomb lattice 
of HgTe nanocrystals (with lattice parameters q = 0.5, a = 5nm) 
compared to the width of the lowest p band, the absolute gap be¬ 
tween the lowest p bands (the gap between band extrema), and the 
vertical gap at T (k = 0 ) between them. £in and Eout are the dielec¬ 
tric constants of the materials composing the lattice (HgTe) and its 
environment, respectively. 


and shape (Supplementary Table 1). 

Under partial filling of the lowest-energy flat p band by 
electrons, fascinating phenomena such as fractional QSHE are 
predicted in presence of strong correlations, when the strength 
of Coulomb interactions between electrons is large compared 
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FIG. 3. Honeycomb lattice of HgTe spheres and its conduction band strncture. a,b, top view of an assembly of spheres connected by 
cylinders, forming a honeycomb lattice of HgTe. Hg atoms are in yellow, Te atoms are in grey. The lattice spacing a, i.e., the center-to-center 
distance between neighbor spheres, is equal to the diameter D of the spheres and d is the diameter of the cylinders. c,d,e, conduction band 
dispersions resulting from the atomistic tight-binding calculation for D = 5.0 nm and djD = 0.3 (c), djD = 0.6 (d), or d/D — 0.8 (e). 
Non-trivial gaps are indicated by pink shaded regions, f, evolution of the gap between the lowest p bands versus d/D for three values of D 
(green triangles: D — 5.9 nm; blue circles: D = 5.0 nm; red crosses: D = 4.1 nm). 


to the band width and smaller or comparable to the energy gap 
between the flat band and the next higher one (for a recent re¬ 
view, see Ref.[3Tl). In addition, it has been shown that also for 
large Coulomb interactions, a fractional Chern insulator phase 
cannot be ruled ouP^ We have thus calculated the on-site 
(t/) and NN (V) electron-electron interaction energies for the 
honeycomb lattice of Fig. [T^, assuming two different dielec¬ 
tric environments around the nanocrystals (Methods). Table[n| 
shows that Coulomb energies are larger than the band width 
and are comparable to the gap between p bands. Therefore, 
the gap sizes open the possibility to the experimental observa¬ 
tion of strongly correlat ed qu antum phases including the long 
sought fractional 

Band structure for other types of HgTe lattices. In order 
to understand the effects of the electronic coupling between 
honeycomb lattice sites, we have studied a second type of 
superlattices with a simpler geometry, consisting of tangent 
spheres connected by horizontal cylinders (Fig. [^,b). The 
HgTe spheres have their [111] axis orthogonal to the lattice 
plane and the cylinders are oriented along (110) directions 
perpendicular to [111]. Figures [^,d,e depict the evolution 
of the band structure with the ratio between the diameters of 
cylinders (d) and spheres (D) (see also Supplementary Fig. 3). 


An increase of d/D induces larger NN hopping terms, broader 
s and p bands, and stronger sp hybridization, as shown by the 
deformation of the s band. It also results in a larger NNN 
term AfgQ, explaining why the non-trivial gap in the s sector 
only exists for d/D > 0.3. On the contrary, topological gaps 
are always present in the p sector, even for small values of 
d/D, because they are mainly determined by the on-site term 
AjgQ. At high values of d/D (> 0.7), the spin splitting of the 
bands in the entire Brillouin zone becomes particularly impor¬ 
tant due to increased Rashba couplings. In general, the lowest 
p band is rather flat and has a maximum separation from the 
next higher one for d/D close to 0.4-0.5. In that case, the 
energy separation can be as large as 35 meV (Fig. H). HgTe 
honeycomb structures with moderate coupling between lattice 
sites should provide the most suitable gap and band widths to 
observe the strongly correlated phases associated to the flat¬ 
ness of the bands. 

In Supplementary Fig. 4, we show further results for a third 
type of honeycomb structure made of overlapping HgTe cylin¬ 
ders parallel to each other. Once again, the band structure can 
be described by the effective model (Supplementary Note 3 
and Supplementary Fig. 5). In that case, the NN couplings are 
even stronger, there is no gap in the s sector due to a large 
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Rashba term, but tbe non-trivial gap above tbe lowest p band 
remains. We can conclude that tbe topological effects in tbe 
p sector are robust against changes in tbe electronic coupling 
and honeycomb period. 


Discussion 

In summary, we have performed atomistic TB calculations 
of the band structure of 2D honeycomb lattices of HgTe. We 
demonstrate that the strong SOC of HgTe combined with the 
honeycomb structure results in several topological phases. 
The calculated band structure can be described by a honey¬ 
comb lattice model as in graphene but including multi-orbital 
degrees of freedom that generate in particular a topologically 
non-trivial flat band. By taking advantage of these features, 
we show that in the same structure not only the QSHE, but po¬ 
tentially also the elusive fractional QSHE could be observed, 
just by varying the electron density. Both topological effects 
turn out to be protected by a gap as large as 35 meV, and can 
thus be observed at high temperatures. Honeycomb superlat¬ 
tices of HgTe are therefore platforms of high interest to study 
electrons on a multi-orbital honeycomb lattice under strong 
SOC. Such structures could be fabricated by nano crysta l self- 
assembly (in a similar way as for PbSe and CdSe^^G^ or by 
a combination of gas-phase deposition and lithography. Our 
results open the path towards high-temperature quantum spin- 
tronics in artificial graphene. 


Methods 

Atomistic TB methodology. The electronic structure of HgTe superlat¬ 
tices is calculated within the TB approximation, as described in detail in 
Ref. 1381 The TB Hamiltonian matrix is written in a basis of atomic orbitals 
{sp^d^s* for each spin orientation) as function of parameters that have been 
obtained by fitting to two reference band structures: Close to the Fermi level, 
we use the k • p band structure of Ref. \M whereas elsewhere we use the 
band structure of Ref, [^ obtained using a quasi-particle self-consistent GW 
approximation in a hybrid scheme. In the present work, we have used the TB 
parameters that give the band structure of HgTe at 300 K. The surfaces of the 
superlattices ai‘e saturated by pseudo-hydrogen atoms that push surface states 
far from the energy regions of interest in this study. Therefore surface states 


do not interfere with edge states predicted in ribbons. Due to the large size 
of the systems that we have studied (up to 5^ 10^ atoms per unit cell), only 
near-gap eigenstates are calculated using the numerical methods described in 
Ref.|^ 

Calculation of the Z 2 topological invariant. The Z 2 topological invari¬ 
ant (spin Chern number Cs) for the bands of interest is calculated using the 
model Hamiltonian following the methodology proposed in Ref. 147 1 and de¬ 
rived from Ref.l4^ This approach works even for systems without inversion 
symmetry, which is the case here. Cg is given by a sum of terms calculated 
on a regular lattice in the Brillouin zone. We have checked that the results 
converge for a mesh denser than 21 x 21 /c vectors. In all cases, the invari¬ 
ants that we have computed for the bands are consistent with the number of 
edge states we observe in the bulk gaps. 

Coulomb interactions. The Coulomb repulsion between electrons in hon¬ 
eycomb lattices of HgTe nanocrystals can be characterized as follows. For 
simplicity, we consider electrons in the s band, since in the case of individ¬ 
ual nanocrystals, it was shown theoreticalljEH and experimentall}^^ that the 
Coulomb integrals ai‘e almost identical for states with s and p envelope func¬ 
tions. The Coulomb interaction associated to electrons on nanocrystals i and 
j is calculated as 

J |VI-i(r)|2|VI,^-(r')|2V(r,r')drdr' (4) 

where T'i(r) is the s state on the nanocrystal i, and V(r, r') is the screened 
Coulomb energy of two interacting electrons at r and r'. Taking into account 
that there is just one (spin-degenerate) s state per nanocrystal, and 
ai‘e simply defined as the components on nanocrystals i and j of the wave- 
function calculated for the lowest s band at F (or another k vector), nor¬ 
malized on their respective nanocrystals. Coulomb matrix elements are de¬ 
composed in the basis of the atomic orbitals and are calculated following the 
methodology described in Ref. im by using usual approximations in the TB 
description, i.e., neglecting overlaps between atomic orbitals and consider¬ 
ing atomic charges as point-like charges. For the potential V, we consider 
two configurations: First, a dielectrically homogeneous system, for which 
V(r, r') = e^/ (£in|r — r'|), where Sin = 14 is the dielectric constant of 
FlgTe; Second, a dielectrically inhomogeneous system, for which V is calcu¬ 
lated by solving the Poisson equation, with the dielectric constant inside (out¬ 
side) the lattice equal to ejn (Sout). We have chosen Sout = 6, a typical value 
taken to simulate the complex dielectric environment around semiconductor 
nanocrystalj^. The on-site (U) and NN {V) terms are presented in Table |n| 
As expected, larger values are obtained for £out = 6 than for £out = £in ■ 
Longer-range Coulomb terms are expected to decay as the inverse of the dis¬ 
tance between nanocrystals. Flowever, it is important to note that these values 
do not take into account the extra screening induced by the electrons filling 
the bands. This could be computed, for example, in the random-phase ap¬ 
proximation, but this is clearly beyond the scope of the present work. Long- 
range interactions will be strongly screened, while short-range ones will be 
only slightly reducecCS In this context, for a band filling of the order of 1/3, 
correlations will be mainly governed by short-range effects. 
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Supplementary Figure 1: Conduction bands in lattices of HgTe nanocrystals. Band dispersion resulting from the atomistic tight-binding 
calculations. Topologically non-trivial gaps are indicated by pink shaded regions. Each nanocrystal has a truncated nanocube shape. 
Truncation factor: q = 0.25 (a,c) or q = 0.45 (b,d). Honeycomb lattice spacing: a = 5.9 nm (a,b) or a = 6.8 nm (c,d). 
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Supplementary Figure 2: Ribbon with broken inversion symmetry. Atomistic tight-binding calculations for a honeycomb ribbon of HgTe 
nanocrystals, with broken inversion symmetry. The nanocrystals have a truncated nanocube shape (truncation factor q — 0.5, honeycomb 
lattice spacing a = 5.0 nm). A single plane of atoms is removed from the right side of the ribbon in order to break the inversion symmetry 
(not shown), f, schematic view of the ribbon. The unit cell of 16 nanocrystals, shown in magenta, is reproduced periodically along the 
direction indicated by the arrow, a, dispersion of the p bands. The position of the bulk bands is indicated by pink vertical bars along the left 
axis. b,c,d,e, 2D plots of the wave functions of four states calculated at fc = 0.3 x 2tt/1 where I is the length of the unit cell. The labels 1 (b), 
2 (c), 3 (d) and 4 (e) refer to the states indicated in a. The plots are restricted to the unit cell of the ribbon. The white dots indicate the atoms. 
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Supplementary Figure 3: Band structure for lattices of spherical nanocrystals. Conduction band dispersion resulting from the atomistic 
tight-binding calculation for honeycomb superlattices of spheres (diameter D — 5 nm) connected by cylinders (diameter d = 0.4_D). a: bulk 
(gap between the lowest p bands = 25 meV). b: armchair nanoribbon composed of sixteen nanocrystals per unit cell. 




Supplementary Figure 4: Band structures for lattices of cylinders, a, top view of an assembly of vertical cylinders forming a honeycomb 
lattice of HgTe (cylinder diameter d — 1.02a = 7.0 nm, layer thickness t = 5.3 nm, a is the lattice spacing). Crystallographic axes of HgTe 
are indicated, b, corresponding conduction band dispersion resulting from the atomistic tight-binding calculation. Non-trivial gaps are 
indicated by pink shaded regions, c, band structure of a zigzag ribbon formed of 12 cylinders per unit cell. 
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Supplementary Figure 5: Effective model for a lattice of cylinders, a, conduction band dispersion calculated using the model 
Hamiltonian, b, band structure of a zigzag ribbon formed by 12 sites per unit cell. 
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Supplementary Table 1: Energy gap in the p bands. Energy gap (meV) between the lowest p bands in honeycomb lattices of HgTe 
nanocrystals for different values of the truncation factor q and the honeycomb lattice spacing a. The same gap predicted for two successive 
values of q means that the honeycomb lattices are identical, the number of atoms per nanocrystal varying discontinuously with q. 


on-site NN hopping Rashba SOC intrinsic SOC 

Es = 266 14.,, = -17.2 = 5.11 Afso = 0.3 

Ep^ = 493 Vppa = 28.9 7pp„ = 4.77 Af’go = 14.1 

Epy = 493 VppT, = —0.6 7pp,r = 0.00 

Ep^ = 698 VsptT = 24.2 

Supplementary Table 2: Parameters of the effective model. Parameters (meV) used for the model for the HgTe superlattice described in 
Supplementary Fig. 5a. Es, Ep^, Ep^, and Ep^ are the onsite energies on the s, px,Py, andpz orbitals, respectively. 14scr, Vppa, Ypp-n-, and 
14po- are the hopping parameters, following the notations of Ref. [1]. 'yssa , 7ppa-, and 7pp,r are the terms describing the Rashba SOC, 
following the same notations. The intrinsic SOC is defined by Afgo and AfgQ on s and p orbitals, respectively. 


Supplementary Note 1 

Influence of size on lattices of HgTe nanocrystals. 

When we vary the lattice spacing a (3-8 nm) and the nanocube truncation {q between 0.25 and 0.5) in lattices of HgTe 
nanocrystals, all band structures have similar behaviour. Examples of band structure are shown in Supplementary Fig. 1. The 
presence of well-separated s and p bands with helical gaps is quite general, even if the width of the gaps may vary substantially 
depending on the geometry of the superlattices. The overall shape of the s bands is always the same, whereas for the p bands 
the variations are more important because not only the nearest- neighbour hopping, but also the respective positions of the p^, 
Py and p^ states depend on nanocrystal size and truncation. However, the lowest p band is always detached from the next higher 
one. The width of the lowest gap in the p bands is given in Supplementary Table 1 for different values of a and q. 
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Supplementary Note 2 

Edge localization in ribbons of HgTe nanocrystal superlattices 

In this section, we present additional results on ribbons made from lattices of HgTe nanocrystals. We consider the same 
nanocrystals as in Figs. 1 and 2 (q = 0.5, a = 5.0 nm) but we investigate a ribbon in which the two edges are not symmetric by 
inversion. First, we discuss the effect of this asymmetry on the band structure. Second, we present plots of the wave functions. 
When the two edges of the ribbon are equivalent by inversion symmetry, the topological edge states on the opposite sides of the 
ribbon are quasi-degenerate for each value of k, the coupling between opposite edge states being negligible (Fig. 2). In order 
to observe the effect of geometry on the results, we have also considered a ribbon in which the inversion symmetry has been 
broken. For that purpose, we have removed all atoms of the last atomic plane on the right side of the ribbon and we have saturated 
the broken bonds with pseudo-hydrogen atoms. Supplementary Fig. 2a shows that the edge states are preserved thanks to their 
topological protection but their degeneracy at a given k is lifted due to the asymmetry between the two sides of the ribbon. The 
2D plots of the wavefunctions of the four states denoted 1^ in Supplementary Fig. 2a for the asymmetric ribbon are shown in 
Supplementary Figs. 2b-e. These states calculated at k = 0.3 x 27r// are strongly localized on the edges of the ribbon. State 3 
is more delocalized than the other three states because it lies very close to the bulk band edge. 


Supplementary Note 3 

Band structures of honeycomb lattices of vertical cylinders 

We have investigated a third type of honeycomb structure composed of HgTe cylinders. The axes of the cylinders are parallel 
to each other and are organized on a honeycomb lattice (Supplementary Fig. 4a). Such structures could be fabricated from 
a HgTe layer, grown for example by gas-phase approaches. The honeycomb nanogeometry is then defined using nanoscale 
lithography. Quite similar band structures are obtained for these lattices (Supplementary Fig. 4b). Once again, s and p-like 
bands can be easily identified. The Rashba SOC is in general much larger than for lattices of nanocrystals due to a stronger 
coupling between neighbouring sites. Interestingly, a very similar behaviour was found for superlattices of spheres connected 
by cylinders when the coupling between neighbouring spheres is strong, i.e., for large values of d/D (Fig. 3). As a consequence 
of the large Rashba SOC, the gap in the s sector is closed. Non- trivial gaps remain in the p bands, for many configurations that 
we have investigated, in spite of larger spin splitting. However, the larger Rashba SOC tends to increase the dispersion of the 
lowest p band. Once again, the results of the atomistic TB calculations are well described by the effective model (Supplementary 
Table 2 and Supplementary Fig. 5a). Only the highest bands of Supplementary Fig. 4b are not reproduced by the effective model 
because they involve higher-energy orbitals which are not considered in the model. The topological properties of the bands are 
demonstrated by edge-state analyses in ribbons, using the atomistic TB calculations (Supplementary Fig. 4c) or the effective TB 
model (Supplementary Fig. 5b) which give results in excellent agreement. The non-trivial topology of the bands is confirmed by 
calculations of the Z 2 topological invariants using the effective-model Hamiltonian. 
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